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FOREWORD 


This  report  was  prepared  by  the  Cornell  Aeronautical  Laboratory,  Incos 
Buffalo,  New  York  as  Com®!!  Aeronautic &1  Laboratory  Report  No.  TB-767-F-7,  under 
USAF  Contrast  No,  AF33(038) -20659 e The  oontraot  was  initiated  under  the  research 
and  development  project,  identified  by  Researoh  and  Development  Order  Number  R-461-1. 
It  was  administered  under  the  direction  of  the  Flight  Researoh  Laboratory,  Wright 
Air  Development  Center  with  Capt.  P,F.  Cerusai  aoting  as  project  engineer. 
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ABSTRACT 


Variations  in  the  longitudinal  stability  of  tne  F-94  airplane  can  be 
achieved  by  automatic  actuation  of  l)  the  elevator  proportional  to  cx-  and  <S$ 
signals  to  alter  the  short  period  damping  and  frequency^?.)  the  pilot's  stick  pro- 
portional to  SF/q  QC  to  change  the  stick  force  and  position  gradients. 


*>)  a 0.43  ft  canard  pitching  control  surface  driven  at  2 deg/sec.  by  u/  and 

f *1 


q 


signals  to  modify  the  phugoid  damping  and  period. 


The  operational  range  of  the  control  equipment  necessary  to  accomplish 
these  variations  in  flight  are  determined.  A specific  sequence  of  calculation  is 
listed  for  obtaining  the  gearing  required  to  obtain  a set  of  flying  qualities. 
Transient  responses  determined  by  analog  computation  are  presented  along  with  a 
phase  diagram  representation  of  the  phugoid  and  short  period  motions. 


PUBLICATION  REVJE.V 


This  report  has  been  reviewed  and  is  approved. 


FOR  THE  COMMANDING  GENERAL? 
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DEFINITION  OF  SYMBOLS 


^ D^-rru)  C*JL 


cL  ^VV\/d  U 
c L^Vn  /ci  q 

cL  %€*/  JL  oc 

1 Se.A6<x 

£•  /d  <x 
DC  ) 
D* 

D o< 
f 

Fs 

' * 

H 

0 

U'B 


Wing  menu  aerodynamic  chord  (ft„). 

Elevator  MAC  (ft.). 

Lift,  drag,  pitching  moment  and  hinge  moment 
coefficients,  respectively. 

Thrust  coefficient,  ■ T 

r 2 qS 

Rate  of  change  of  pitching  moment  -with  non- 
dimensional  airspeed,  dC^d  (u/V). 


Rate  of  ohange  of  pitohing  moment  with  non- 
dimensional  rate  of  change  of  airspeed. 


dCm/d 


Rate  of  ohange  of  pitohing  moment  with  airspeed  (rad/mph) 

Rate  of  ohange  of  pitohing  moment  with  rate  of  ohange 
of  airspeed,  dCj^^dvy/d^rad/mpiyBeo) 

Rate  of  ohange  of  elevator  angle  with  angle  of  attack 

Rate  of  ohange  of  elevator  angle  with  rate  of  change 
of  angle  of  attack,  (deg/deg/Beo) 

Rate  of  ohange  of  downwash  angle  with  angle  of  attack 


d ( )/d(t/tr),  D2  - d2  ( )/d  {t/t,)Z 


Rate  of  change  of  non-dimensional  airspeed  with  non-dimensional 
time,  d(u/v/d(t/tj) 


Rate  of  ohange  of  angle  of  attack  with  non-dimensional 
time,  doc/d  (t/t?) 


Short  period  frequenoy  (ops) 

Stick  force  (its.) 

Acceleration  of  gravity 
Hinge  women t (ft.  lbs.) 
Non-dimensional  pitohing  inertia. 
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p 

f 

Vh 

£ 

s* 

•D 

T 

a 

9 

u 

v 

w 


Jf 

Ss 

M 

5P 

S* 

Q 

f 


Tail  length  (ft.) 

Normal  acceleration  increment  in  g’s. 

Phugcid  period  (seo.) 

Dynamic  pressure  (psf.) 

Horizontal  tail  dynamic  pressure  (psf.) 

Wing  area  (ft2) 

Elevator  area  (ft^) 

Time  (seo.) 

Thrust  (lbs.) 

Forward  speed  increment 

Forward  velocity  gust  disturbance,  (mph). 

Rate  of  change  of  airspeed  with  time,  du/dt  (mph/eeo.) 
Forward  speed,  V*  denotes  indicated  airspeed  (mph.) 
Gross  weight  (lbs.) 

Angle  of  attaok  (deg.) 

Rate  of  change  of  <x  with  time,  do</dt  (deg/sec.) 
Auxiliary  surface  deflection  (deg.) 

Elovator  deflection  (deg.) 

Flap  deflection  (deg.) 

Control  stiok  deflection  (deg.) 

Elevator  tab  deflection  (deg.) 

) Increment  in  ( ) . 

phugoid  damping  ratio. 

Short  period  damping  ratio. 

Angle  of  pitch. 

Relative  density, 

Air  density  (lbs/ft®) 

TMrnn  iv«i  -h  „ /j/ '/o  i 9 )/ 

“ ’ * f t n 
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INTRODUCTION 


The  Flight  Research  Laboratory  of  the  Wright  Air  Development  Center  has 
instituted  a program  with  the  Cornell  Aeronautical  Laboratory  to  obtain  actual 
flight  test  data  on  the  optimum  and  minimum  flyable  longitudinal  stability  and 
oontrol  oharaoteristios  for  fighter  and  bomber  airplanesc  This  type  of  information 
has  recently  become  of  great  design  importance  with  the  advent  of  praotioal  servo- 
mechanisms  for  the  addition  of  artificial  stability  to  airplanes?  also,  this 
information  should  be  useful  to  those  oharged  with  the  responsibility  for  estab- 
lishing handling  qualities  specifications. 

Two  airplanes  are  being  used  for  the  evaluations-  one,  a B-26  light 
bomber?  the  seoond  an  F-94  jet  fighter*  The  elevators  of  these  airplanes  are 
driven  by  irreversible  hydraulic  servos  in  response  to  oontrol  signals  supplied 
by  the  pilot  and  signals  provided  by  artificial  stability  pickups*  Ihe  control 
sticks  are  driven  by  a seoond  servo  in  response  to  pilot  applied  control  force,  in 
a maimer  closely  simulating  the  natural  airplane's  control  forces.  A small  auxiliary 
pitching  surface  is  driven  by  an  eleotrio  servo  motor  for  phugoid  oontrol.  By 
adjusting  the  gains  of  the  various  channels  of  this  equipment,  the  following  parama- 
ters  of  longitudinal  stability  and  control  oan  be  varied?  phugoid  mode  period  and 
damping,  short  period  mode  period  and  damping,  static  elevator  to  trim  vs.  and  g, 

and  static  stick  force  vs.  Cl  and  g.  The  extremes  of  stability  and  control  "that 
oan  easily  be  simulated  and  evaluated  could  not  safely  and  economically  be  obtained 
in  any  other  way* 

Reference  (1)  presented  the  theoretical  analysis  on  the  B-26  airplane. 

This  report  contains  the  supplementary  analysis  of  the  artifioial  oontrol  necessary 
to  provide  a wide  range  of  flying  qualities  for  the  F-94  airplane* 
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ANALYSIS 


BASIC  FLYING  QUALITIES 

The  eight  flying  qualities  which  are  considered  to  be  basic  in  a pilot's 
evaluation  of  an  airplane's  longitudinal  motion  were  Hated  in  Inference  (l)  as 
follows* 

1.  Frequenoy  of  the  "short  period"  mode 

2.  Damping  of  the  "short  period"  mode 

S.  Period  of  the  "phugoid"  mode 

4.  D&nping  of  the  "phugoid”  mode 

5.  Control  deflection  to  trim  vb.  speed  for  lg  flight 

6.  Control  defleotion  per  "g"  normal  acceleration 

7.  Control  force  to  trim  vs.  speed  for  lg  flight 

8.  Control  force  per  "g" 

The  subsequent  analysis  will  show  how  these  baBio  quantities  will  be 
varied  on  the  F-94  airplane. 


Short  Period 


(Ref.  1) 


The  short  period  oharaoteristios  can  be  determined  frisa  the  formulae* 

£ 


f = 


100 


JL 


^8 


The* 


^4-  a*  ^")*i  a - 

— r-" 


Z s* 


(100) 


it 


2T^( 


_ 

- — ~ 


!"o W 


eps 


(1) 


{%  critical) 


The  normal  stability  derivatives  of  the  F-94  are  listed  in  Table  I.  The  frequency 
and  per  cent  critical  damping  in  the  short  period  mode  are  listed. in  Table  II  for 
five  possible  flight  conditions. 

Possible  artificial  control  of  the  short  period  at  292  mph  indicated 
airspeed  and  20,000  ft.  altitude  is  listed  in  Table  III.  it  is  noted  that  either 


0«  2^  ohanges  the  effec- 


ACm-.  «r  Cm2  o*a  be  used  to  vary  the  frequency. 

* jD  Q 

tive  pitching  inertia  of  the  airplane,  and  its  use  yields  only  a limited  range  of 
frequenoy  control,  especially  with  large  statio  margins.  ZiCja^  waB  chosen, 

therefore,  bo  allow  wider  short  period  frequenoy  variations. 

The  short  period  damping  can  be  controlled  by  either  A Cjj^  0T  ^cmq 
as  ssen  in  Table  III.  Acm  » however,  affects  the  phugoid  mode  and. the  maneuvering 

stability.  AC™  was  chosen,  therefore,  in  order  to  isolate  the  effects  of  each 
in- 
variable on  the  eight  basic  flying  qualities.  It  should  be  noted  that  either  Ac, 
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or  Ac  ohanges  both  the  frequency  end  the  dumping  of  the  short  P»rlo 

’6“  , f0r  these  ertlflel.1  control,  he.  been  established  as 

Thus,  a sequenoe  of  settings  lor 

follows* 


1.  Choc,  desired  f Bid  °f  ■»"*  S”**’ 

2.  Calculate  dSe/dex  OF^cm^  required  from? 


(Ref.  1) 

sL  -4^2 


r C-usC~*< 

>v" 


3.  Calculate,  dfe/d*  or  Ac^  fromJ 


j£jt 


C. „ C-vy, , 


2 **  "s  ' 
(rad/rad/sec) 


— 


;!"a  <2> 


-E9s 


S'  "*(?> 

c cTfl 


Phugoid 

Th,  phugoid  ohereoterietios  of  the  normal  *-94  «rpl»  can  be  determined 
from  the  formulae i 

<R«f- 11  r . c,s  Cn,  lr 


ysLk- 


*RCL  £*0, 


h = ^ t Cj^r* 

JX  the  static  margin  is  greater  than  10jS.  the  phugoid  psriod  and  doping  are 
closely  approximated  by  the  formulae! 

P - .202  (»•»•) 

100  C - 70.7  W critical) 
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The  normal  phugoid  characteristics  of  the  P-94  should  be  calculated  from 
the  former  set  of  equations  as  the  statio  margin  is  close  to  5.5$  for  the  usual 
flight  conditions.  Table  II  lists  the  period  and  damping  for  five  possible  flight 
configurations. 


Possible  artificial  control  of  the  phugoid  at  292  mph  indicated  airspeed 
and  20,000  ft.  altitude  is  listed  in  Table  IV.  The  purpose  of  this  program  is  to 
provide  a wide  range  of  control  over  the  flying  qualities.  Variations  in  thrust 
(or  C^t,  ) would  not  provide  this  extreme  range  in  controlling  large  amplitude  dis- 
turbances. This  oan  be  seen  for  example,  if  a = -.17  is  assumed  to  add 

^ Cl 

50$  critical  damping  to  tho  phugoid.  Tf  the  airspeed  varied  by  20  mph  ind„  a change 
in  thrust  of  1220  lbs  would  be  needed,  to  achieve  this  damping.  Thus,  it  was  not 
considered  practical  in  this  case  to  utilize  thrust  variations  for  phugoid  control. 

Pitching  moment  derivatives  proportional  to  u or  D^u  oan  bo  used  to 
vary  the  phugoid  period,  » however,  affects  the  short  period  mode,  bo  that 

Cmu  has  been  chosen  for  artificial  control. 


Either  Cn^  or  _ will  control  the  phugoid  damping.  Since  C. 


m 


m. 


would  also  vary  the  static  stability  of  the  aircraft,  C was  chosen.  It  should 

“bu 

be  noted  that  boih  ^ Cjn^  should  be  used  simultaneously  to  achieve  a 

desired  period  and  damping,  sinoe  either  derivative  affects  both  period  and  damping. 


Auxiliary  Surface 

It  is  noted  in  Table  IV  that  a dSe/du  of  -.0059  deg/mph  ind.  would  be 

required  to  reduce  the  phugoid  period  to  62.8  sec.  This  gearing  would  demand  a 
positioning  accuracy  of  0.0069  deg.  elevator  if  the  threshold  of  the  instrumentation 
is  assumed  to  be  1 mph  airspeed  variation.  Therefore,  the  elevator  can  not  be 
used  to  provide  the  low  pitching  moments  required  to  control  the  phugoid.  The 
design  of  a small  pitching  control  auxiliary  surface  is  now  considered. 


The  maximum  pitching  moment  neoessary  will  in  general  determine  the  di- 
mensions of  the  auxiliary  surface.  In  this  case,  however,  the  design  was  determined 
by  calculating  the  minimum  moment  deBired  and  then  multiplying  this  value  by  100  to 
yield  a practioal  gearing  range  for  instrumentation.  The  minimum  Cm  would  be 
required  to  add  10$  critical  damping  to  the  phugoid  at  the  Vi  = 292  mph,  H • 20,000  f t.  3 

dCffl/dCL  = -.055  condition  when  the  disturbance  iB  of  1 mph  amplitude. 


&n!  °"b“DU 


°“!>u  ' TOT  4<7s 

(Ref.  1) 


2dCffl/dCL  (dCJjn/dCL  + 
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034 


-m, 


Du. 


4(.10)  |I(-.055)  (-e066l)j 


Du 


V 


JIEmin. 


fZlT'CL  1 2/T(2.4)  _ .000583 

ITT  757  "W.z-  — 

,034(. 000583)  - .00198  x 10~2 


I?  * 100  to  1 ratio  is  assumed,  for  the  sensitivity  range,  a design  Cm  of  .002 
will  be  determined. 

A defleotion  range  of  1 io  deg.  was  assumed  to  give  a Cm  of  .0002 

Sa. 

l/dego  The  dimensions  of  the  surface  were  determined  from  the  formula* 


For  an  assumed  aspect  ratio  of  2.8?  for  each  side  of  the  8urfa.ce,  C j,  will  be 

3.06.  The  surface  will  be  mounted  in  the  nose  14  ft.  from  the  c.g.-of  the  airplane. 
Thus* 


Sa  * — .0002  (57 o3)  — , o43  ft2 
„ / 14  a 1 


3.06 


(7777) 


737T 


The  dimensions  of  the  aurfaoe  are  then* 


spaa  on  one  side  = 9 3/8  in. 

chord  =3l/4  in. 

The  formulae  for  dSa./du  and  dS«/du  in  terms  of  the  phugoid  damping  and 


period  ares 
d&'K?  ' 

■t* =•  t ■»  * 1 • cm  r— 7— 3 


dt  <4 


where  Q 


wg 


'So, 


r /2  'T'T  I 

L^r  ~ l?TT 


2^7" 


(5a) 


deg/mph  ind. 


= C„  4-  C 

^ ex  I 


cLSe/tlc 


h 


T 


Cyr, 


°<T  i~ 


2/« 
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<i<4  ^ / c°  ^ 5^7, 


AS,  ~h  C*  //-£ess_i 
V'  C,  / 


<L 


r c.. 


? elpr 


(6a) 


where 


>>i 


=c 


_ (r^ui 

£ 


D<v  f — >* 


C* 

*£o<_ 


dSm/dU 


-9, 


ii‘ 

fir  I 
L-  A 


deg/mph 

ind/seo. 


and 


$ = C», 


C±af  y ,-  f"  d 

r^r^^soinr 


The  sequence  of  calculations  is  important.  Steps  (1),  (2)  and  (3) 
should  be  completed  first  for  the  short  period  analysis  and  then; 

4.  Choose  ^ «d  P of  the  phugoid  desired. 

5.  Calculate  dSeo/du  or  as  above. 

6.  Calculate  dSau/du  or  C_  as  above. 


It  was  noted  in  Reference  (l)  that  the  signals  A u/q  would  be 

provided  instead  of  u and  u in  order  to  vary  the  phugoid  in  a "natural"  way, 
i.e.  independent  of  speed.  The  signal  A -Jc,  was  provided  by  a dividing  network  as 
described  in  Reference  (2)  whioh  allowed  the  full  range  of  q while  discerning  a 
1 mph  variation.  The  q range  of  the  F-94  is  greater  than  that  of  the  B-26  and 
with  the  use  of  a larger  range  gage,  the  discrimination  can  no  longer  be  held  near 
1 mph  (or  1 paf.). 


A Aq/q  signal  will  be  used  for  the  F-94  instead.  This  signal  can  be 
measured  accurately  by  differential  pressure  gages  on  total  and  statio  lines  and 
with  a dividing  network  sb  before.  The  relations  between  A -A  and  A q/q  can  be 
noted  as*  » 


H-p  A H-  hH-*-  t 

Thus  it  is  apparent  that  Aq/q  is  proportional  to  <4  (/f)  and  that  the  Aq/q  signal 
will  be  inversely  proportional  to  the  equilibrium  dynamic  pressure.  This  signal 
must  then  be  adjusted  for  the  trim  speed  of  the  test. 


The  auxiliary  surface  will  then  be  controlled  by  signals  from Aq/q  and 
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8a.  Caloulat-e  d5e/d5s  (If  d§fi-/dVi  is  chosen) from 

44^  = 

d*  V^6  £)** ok  s6  V^/ 

or  8b.  Calculate  d$e/dS5  (if  d§e./dn  is  ohoBen)  from* 

d'Ss  c^e 


9.  Choose  d%/d?i  and  d^/dn  desired 
10.  Calculate  d^^/doc  from 


deg/deg  (8a) 


deg/deg  (8b) 


sLSjSu-.-  _SjU&fe- L !-  ifel  - ~ j hr  4^- 

<*'  <*  c, 

Sf  5$ 


^Es/d'iu  //■ 


(at  trim  speed  in  Ig  flight) 
11.  Calculate  dFs/q/dS>  from 


deg/deg 


JbFs! 

cb%s 

- #- 

igL-  Se* 

dfFsfp 

r l£k_ 

c^s 

c l Ss 

e)  VJo 

2 WJ  s 

'$S 


■« 


1 


VKC 


Ss 

(ft?/ deg) 


(11) 


PREDICTING  FLYING  QUALITIES 

The  steps  to  predict  the  flying  qualities  are  listed  below  as  a summary 
of  the  preoeding  analysis i 

1.  Choose  jP  and  Cjs  of  short  period  desired. 

2.  Find  dS«,/d*i  from  Equation  (2). 

3.  With  this  d«e/do<)  use  Equation  (3)  to  determine  d§*/de<  or  A Cjg^. 

4.  Choose  P and  of  phugoid  desired. 

5.  With  dSft/do,  and  di’e/de?  found  above,  use  Equation  (6a)  and  (5b) 

to  determine  or  Cm^, 

6.  With  d #«.  / do* j dSs  /da^  and  dC«,/du  determined  use  Equation  (6a)  and 
(6b)  to  find  d&ftz/du/q  or  _ 

- 7 - 
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7.  Choose  either  dSc/dY±  or  d5®/dn  « 

0*.  Witii  ^Ca#<  and  Cn^  known,  use  Equation  (8a)  to  find  Cmgs  or  dSe/d§s 
if  dSe/dVj,  is  ohoaeno 

or  8b.  With  and  Ca^  known,  uaa  Equation  (8b)  to  find 

if  dSfi/dn  is  oho san. 

9o  Choose  dFs/dVi  and  dFg/dn  desired. 

10.  With  ^Cma^  known,  use  Equation  (10)  to  determine  dSs /de^  or 

<i-»AhV 

11.  With  and  /d&p  known,  us©  Equation  (ll)  to  determine 

0hs*  or£f^‘ 

KAMGK  OF  VALUES 

The  gearings  dSe/dof)  d$e/do(  d’&'u/d^j^-j  dS^/d^j 
dS s /do<  and  dFg/c^/d Ss  have  been  chosen  for  automatic  controls 
values  of  these  parameters  will  next  be  determined. 

Elevator  Servo 

Modifications  to  the  normal  F-94  airplane  will  be  made  to  aotuate  the 
elevator  with  no  meohanioal  connection  to  the  pilot5  a sticks  It  should  be  noted 
here  that  a stiok  will  be  added  in  the  F-94  aft  cockpit  to  allow  the  stand-by 
pilot  to  maintain  direct  mechanical  oontrol  of  the  elevator  surface. 

An  electronic-hydraulic  servo  system  similar  to  that  described  in 
Reference  (2:)  will  be  used  to  aotuate  the  elevator  surface.  Full  elevator  travel 
of  54  deg.  will  be  provided  with  6 in.  servo  linear  travelo  The  design  resolution 
of  the  servo  is  l/6  of  1%  or  .01  in  whioh  is  equivalent  to  .09  deg.  The  signals 
to  the  elevator  servo  are  shown  be  lows 


dfa/dSs  > 

The  ranges  of 


The  equipment  will  be  designed  for  a short  period  frequency  extremes  of  0.1  to 
1.0  cps  and  for  an  angle  of  attack  of  from  ±0.07  to  ±io  deg. 


change. 


The  signal  d%e/d c*  will  be  uBed  to  provide  a minimum  1 % static  margin 
Thus,  from  Equation  (2)  % 


c 


vw 


wZ 


+ 

Cj 


5.9  (.01)  ■ 
.945 


-.0625  deg/deg 
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“•  0 CO 


. .v_  neriod  be  doubled  at  the  high  atatio 

•S’SSL.  «•»»  °f  45e/d-  - 

3 V 

cr 


dS^ 

a c*. 


“3.  - a(b.9)  ( -•11)  . 2,06  deg/deg 
" “945 


5e 


In  order  to  halve  the  frequency  at  this  conditions 

_dSe.  _ —5-  m “eB1  d#s/d*B 

- , .1  a will  provide  « mini™-  U*  orltlo.1  demping  ohenS. 

*“  .t tl  - mph'lnd.  end  2)2  .tail.  m«gln. 

in  the  short  period  at  Vi  f 

From  Equation  (3)» 


_ _ Z^A^f 

<Z'  o<  v' 


/ _ 

( 3f,j. 


%*■) 


■4- 


. + 2(2,40).10  fiiglrR.flU.02  + .0111)]  - - «026  deg/deg/eec. 

. _ ' to  incre.no  the  .tort'  dcmpln.  to  1^^ 

r“  aw  % i»a-  “d  io!S  ■t*tic  “"s1”-  ’ 

and  d Je  /d  of  lo* 

<££e  2t2.«)i;°M 

In  order  to  bring  the  damping  to  aero* 


|i|i  (5.9)(.10  +.0111) 


2 - .48  deg/ deg/ seo 


Jgfr  - M||*l  (..496)  [.09s]  - -.25  d.g/d.g/..o 

^ , (O-  ,, _d.d_V  C M 


r«-y-  545  ' ' j 

stiok  noeition  gradient  with  airspeed  and 
The  signal  S ^ d §e/dVi  and  djfi/dn  from  slightly 

aias.'B..-.  •y^Siiartusr 

Minimi—  - t 0-1 

■ - d°  dse/a=S 

Stick  Servo 

. i -tiok  will  also  be  actuated  by  an  eleotronio-hyd  aulio 
“ £3  dog.  will  t.  prorided  'ey  6 in.  -«*- 

SSlfS’.iSSi.  to  the  .tide  ..rwo  .r.  oho»  below. 

F*/$> 


a<_ 

St 


— — - 

c a 

STICK 

©>S 

SERVO 

SYSTEM 

9 - 
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The  stiok  foroe  gradients  will  be  varied  by  signals  from  F<t/q  and  eX,  while  -the 

St  signal  will  be  used  to  trim  the  stick  to  zero  foroe . 'The  strain  gages  on  the 
aluminum  stiok  will  measure  strain  equivalent  to  10,000  pci  stress  when,  the  etiok 
foroe  is  100  lbs.  The  dynamic  pressure  range  for  the  F*=94  will  be  from  46  to  340  psf. 


The  stick  foroe  per  speed  change  (dFB/dVi)  Qaa  k®  determined  from 
Equation  (7)  asi 

Cv 


dF, 

dVi 


dCa 


- 0, 


%t 


dCL  Ch 


8 


5he  atiok  foroe  per  stiok  defleotion  oaa  be  expressed  as  a total  derivative  in 
terms  of  the  servo  gains  ass 

® ■&  ■ *©) « (2r)  * (£)  «(£) 


Total 


servo 


servo  airplane  servo 


The  servo  gain  dFg/dg  8 is  then  proportional  to  the  total  derivative  and  is  equiva- 
lent from  Equation  (li)  toj 

d**8  Cu 

Ch«» 

Variations  in  ChgB  or  dFa/dfis  will  provide  control  over  the  flying  quality  dF8/dVio 

If  a d?e/dVi  range  from  slightly  unstable  to  five  times  the  normal  .stable  value  is 
required  and  a faotor  of  two  is  allowed  for  d$e/d$B  or  dgs/do<  variations,  the  .range 

of  0*.  and  dFB/dge  will  be* 


Maximum  Gh^g  “ 10  °hie 


-2 


A practical  minimum  range  would  be  8 
ch£B  m *1  Chje 


or  dFg/dgg  ■ 107  lbs/deg. 

-21  lbs/deg. 


or  dFB/d£B  *»  ±1,1  lbs/deg« 


+5 


For  a Ch^s  range  of  t, 2 C^r,  to  -2  ch£e  afc  the  extrem*  limits  of  q,  the 


normal 


maximum  dFg/dSs  will  be 

dFs/dgg  ■ 11  lbs/deg  at  q “ 46  psf 
dFs/dgg  “ 84  lbs/deg  at  q - 340  psf 
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The  signal  d«/d£s  will  be  used  to  provide  the  dFg/dn  gradients  required. 


At  q » 218  psf  the  range  of  stick  foroe  per 
minimum  amount  of  control  was  calculated  as 


g will  be  from  -20  to  +50  lbs/g. 
a der/dgs  equivalent  to  a C}^  of 


The 


t.10  l/rad.  This  is  equivalent  to* 


dS's 

def 


i.020  deg/deg 


The  stiok  foroe  per  g range  can  be  provided  if  be  -1.06  and  +1.30  l/rad. 

respectively  for  dFs/dn  of  -20  and  +50  Ibs/g.  Including  a factor  of  five  as  above 
for  Ch£g  variations,  dgg/d^  will  be* 


Maximum  negative 


foe 

d«< 

Maximum  positive  = 


1.30  (.475) 

•482  W) 

L 1.05  (.475) 

^ Wf) 


a -6.4  deg/deg 
+5.2  deg/deg 


Auxiliary  Surface 

The  auxiliary  surface  will  be  actuated  by  a sm&\ll  electric  servo  motor 
geared  down  for  a maximum  2 deg/seo  surface  rate.  HiIb  velocity  limit  will  allow 
the  surface  to  control  the  phugoid  without  causing  undesirable  normal  acceleration 
response  to  u gust  inputs  as  shown  later.  This  Blow  rate  will  be  sufficient  for 
phugoid  control  ub  a full  range  of  +10  to  -10  deg.  surface  travel  will  take  20  sec. 
or  about  one  quarter  of  the  phugoid  period  at  292  mph  ind.  The  signals  to  the 
surface  are  hewn  below* 


AI3C.  SURFACE 
SERVO 
SYSTEM 


s*. 


The  auxiliary  surface  position  signal  proportional  toAq/q  will  be  pro- 
vided to  vary  the  phugoid  period.  The  surface  will  have  a t10  degree  position  range. 
The  largest  control  range  was  assumed  at  the  Vi  **  292  mph  condition.  A minimum 
A<1  of  1 psf  (or  1 mph)  and  a maximum  Aq  of  129  psf  were  assumed.  At  the  extreme 
q range,  aAq  of  -5  psf  at  q » 46  psf  and  a Aq  of  ±24  pBf  at  q - 340  psf  were 
assumed. 


The  minimum  signal  will  provide  for  Yfc  ohange  in  the  equiv.  static 

nargin.  Thus  from  Equations  (5a)  and  (5b) j 

dga  = Gmu 
d ^ 2fimga 


2£El_  (a  iSsU  + -264  , 01)  B + 

2CmSa  l dCj  .0002  U01'  _13« 


2 deg< 
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Thu  maximum  d£*/d^r  signal  will  allow  a halving  in  the  phugoid  period 
at  6,695  static  margin  and  Vj.  * 292  mph« 


!l£  “ - ^64^  (+.166)  - 218  deg. 


2Cffl<a 


The  maximum  negative  value  is  found  if  the  period  is  to  be  doubled. 


<4a  . 


^ji|)(-.055)  - -S5  an 


The  signal  proportional  to  u/q  will  be  provided  to  change  the  phugoid 
damping  a minimum  of  5$  oritioal  at  5.5/5  statio  margin  and  7^  * 292  mph.  From 
Equations  (6a)  and  (6b)  - 

+ 


d«Sa  q t 


i 


*£»,  in) 

dSL  TJa) 


d^m  / 

dClidCL 


f.  ” v.  j_. - - ' J J 

■ i§  [*(—)(-««)  7 ■ n. 


„„  deg  sec  psi 

06 mS& iar.'” 


The  maximum  sensitivity  of  d§*/d-£-  139  required  to  increase  the  damping  to 

70^5  oritioal.  Bauss 


dga  _ 218  (2.4 


)4(.70  - .059) 

V ^ L 

In  order  to  reduce  the  damping  to  -2 ($  oritioal 8 


j~2( -.055)  (-.0661)  ^ 


13.6  dog  psi 

mrih~ 


dSa  _ 218  (2.4)  ^/  or^ 

144  (2W^SM  J4(-.269) 


2 (=.055) (-.0661) 


-5.7  deg  aeo  psi 
mph.  ind. 
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Summary 


Derivative 

Elevator  Servo 

4 ,0625  to 

+ 2.06  deg/deg, 

- .51 

+ ,025  to 

+ 0i*8  deg/deg/sec. 

- .23 

/^S 

1 0.1  to 

1 10  deg/deg. 

Stick  Servo 

t//=s/<a/Ss 

* 1.1  to 

4 10?  lbs/deg. 
- 21 

t .02  to 

4 5,2 

- 6,li  deg/deg. 

Auxiliary  Surface 
Servo 

t 13.2  to 

4 218 

- 55  deg. 

t 1.06  to 

4 13.67deg/sec/psi  \ 
- 5.7'  mph  ind.  / 
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TRANSIENT  RESPONSE 


The  response  of  the  F-9U  at  292  mph  indicated  airspeed  and  20,000  ft. 
altitude  was  determined  by  analog  computation  for  three  disturbance  conditions.  A 
step  Cm  was  introduced  on  the  analog  by  either  a sudden  application  of  1 deg. 
elevator  or  10  dag.  auxiliary  surface  deflection.  Ousts  in  and e<  were  entered 
as  initial  conditions.  A description  of  the  computational  approach  can  be  found 
in  the  appendix* 


Pitching  Moment  Steps 

Figure  (1)  shows  the  F=9U  response  to  a 1 deg.  siep  elevator  deflection 
for  the  normal  condition  of  \ ■ 292  mph,  20,000  ft.  altitude,  5,5%  static  margin 
and  gross  weight  of  13,611;  lbs.  The  heavily  damped  short  period  is  evident  in  the 
slight  anglo  of  attack  overshoot.  The  low  phugoid  damping  is  also  apparent. 

The  effects  of  Zl  and  Z^  cm^  are  seen  in  Figure  (2).  The 

damping  of  tho  short  period  is  increased  with  negative  '»  The  elevator 

trace  shows  that  about  0.2  deg  of  elevator  control  would  be  needed  with  this 
artificial  control  for  1 deg  step  elevator  input.  If  negative  Zs.  CL.  . is 

added  the  spring  constant  of  the  motion  changes  with  a corresponding  increase  in 
frequency.  The  1 deg.  elevator  step  now  creates  less  and  response  due  to 

the  increased  stability.  Over  0.7  deg.  elevator  control  is  required  with  this 
artificial  derivative. 


Figure  (3)  reveals  the  phugoid  mode  variations  with  C 


and  C™ 

U£C* 


Positive  Gju^^ 


of  .171  adds  $0%  critical  damping  and  requires  almost  5 deg.  of 


auxiliary  surface  control  for  each  10  deg.  of  excitation.  The  decreased  period 
due  to  C ■ .029  is  evident  in  all  the  traces.  This  artificial  control 


demands  almost  as  wide  a variation  in  as  the  excitation  itself  as  seen  in 

the  auxiliary  surface  trace. 


Forward  Velocity  Gust 


Tho  re  sponge  of  the  normal  F-9li  to  a 10  mph  ind.  forward  velocity  gust 
is  seen  in  Figure  (U).  This  gust  as  represented  by  an  initial  condition  on 
causes  an  initial  value  of  n at  zero  time.  After  the  short  period  ot  peak,  -the 
motion  is  characterized  by  the  lightly  damped  phugoid  mode.  If  .029  Cm,  is 


added  by  artificial  control  the  gust  reponse  changes  to  that  shown  in  Figure  (9), 
The  g response  increases  somewhat  after  its  initial  value.  Approximately  7 deg, 
of  auxiliary  surface  deflection  would  be  required  to  provide  this  artificial 


G, 


m . 


in  response  to  a 10  mph  gust. 
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Figure  (6)  shows  the  effect  of  0 on  the  normal  acceleration 

"D  u. 

response  due  to  a forward  velocity  gust.  In  the  normal  case  at  the  left,  the 

in  itial  g decreases  steadily  after  the  short  period  dip.  If  a (L  of  .171 

D cu 

is  added  to  the  airplane,  the  response  shows  a sharp  peak  about  seven  times  the 
initial  normal  acceleration  amplitude.  This  effect  is  due  to  the  impulse  in 
Cjj  created  by  use  of  » as  explained  in  Reference  (1).  The  use  of  a 

velocity  limiter  on  the  right  hand  graph  of  Figure  (6)  reveals  the  close  to 
normal  response  obtained  by  limiting  the  rate  of  buildup  of  to  2 deg/sec- 0 

Angle  of  Attack  Gust 


The  effects  of  A Cja  . and  A Cm.  on  the  F-9U  response  to  an 

angle  of  attack  gust  is  seen  in  Figure  (7).  The  magnitude  of  the  gust  input 
varies  from  „J$2  deg.  for  the  normal  airplane  and  for  As  CL  , , to  1 deg. 

for  the  A condition.  This  limitation  arose  from  computer  overload 

problems.  Overshoot  in  angle  of  attack,  pitch  rate  and  normal  acceleration  is 
evident  in  all  but  the  A C condition.  Mien  dO%  critical  damping  of 

the  short  period  is  caused  with  A the  aircraft  responds  smoothly 

to  return  to  the  equilibrium  condition.  The  use  of  A C. 


-,0U3  would 


require  about  0.2  deg.  elevator  per  deg.  angle  of  attack  gust,  while  use  of 
A Cm  . * -*389  requires  0.6  deg.  JT  par  deg.  ^ gust.  The  increase 

in  the  short  period  frequency  by  A CIil  . is  evident  in  the  faster  initial 

response  to  the  gust  input.. 


It  should  be  noted  that  the  transient  curves  presented  were  transcribed 
directly  from  Brush  recorder  analog  results.  Thus,  the  time  scale  is  correct 
only  at  the  equilibrium  value.  A similar  time  at  any  amplitude  can  be  found  on 
a circular  arc  of  3 inch  radius.  It  was  not  deemed  important  to  alter  this 
time  scale  to  the  usual  rectangular  coordinates. 


VECTOR  PHASE  DIAGRAMS 


The  homogeneous  equations  of  motion  of  the  aircraft  can  be  solved  for 
the  vector  balance  of  forces  and  inertias  required  to  maintain  equilibrium  at 
any  frequency  of  oscillation.  This  vector  balance  is  illustrated  in  Reference 
(3).  The  method  of  obtaining  the  phasings  of  the  motions  in  the  phugoid 
and  Mia  siiui' v period  is  presented  in  ths  Appendix.. 
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Short  Period 


The  short  period  phasing  for  the  F-9UA  airplane  at  « 2y2  raph, 

20,000  ft.  altitude  and  $.$%  static  margin  is  presented  in  Figure  (8).  It 

can  be  seen  that  for  the  normal  airplane  the  D © response  will  be  7 times 

the  #<  response  and  87 • 2°  ahead  of  ^ > The  amplitude  and  phase  of  D&4  in 

relation  to«a^  indicates  the  ,31U  cps  frequency  and  k9«$%  critical  damping  of 
the  motion.  The  unimportance  of  ct  in  the  short  period  is  apparent. 

The  effects  of  A and  ^ on  the  motions  are  illustrated 

in  Figure  (8) . Both  artificial  derivatives  raise  the  frequency  to  .1*7$  cps, 
as  can  be  noted  in  the  similar  amplitude  of  the  D^,  vectors.  The  use  of 
Cm  will  lower  the  damping,  however,  as  seen  by  the  smaller  lead  angle 

between  D<A  and  * The  D & response  due  to  zit  will  be  almost 

in  phase  with  the  normal  D©  vector,  while  Cm  2 will  cause  D©  to 

lead  <s*  by  101.3°. 


Figure  (9)  shows  the  effects  of  Cm_  and  A C, 


■m  on  the  responses. 

q 


Both  derivatives  increase  the  damping  to  19-$%  critical,  but  A Cm^  also  causes 

an  increase  in  frequency  as  seen  by  the  larger  amplitude  of  the  D<A  vector. 

The  increase  in  B & lead  over  ^ is  apparent  for  both  artificial  derivatives, 
with  A (V,  also  decreasing  the  amplitude  of  the  D O vector  considerably. 
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Phugoid 


The  phugoid  phasing  for  the  F— 9UA  condition  is  shown  in  Figure  (10). 
The  low  damping  (5.9^  critical)  is  apparent  in  the  0 u.  lead  of  u.  by  only 
93. It0.  The  high  period  of  88,6  sec.  is  noted  in  the  length  of  the  Dm  vector. 
The  phugoid  motion  occurs  at  almost  constant  angle  of  att^gk,  as  noted  by  the 
minute  <=»<*  vector,  while  the  pitch  angle  variation  is  quite  important  being 


1.28  times  the  amplitude  of  4 or 


Figure  (10)  also  illustrates  the  effects  of  C > t ^ 


and  lagging  t(  by  95>°. 

and  C 


7^  Du 


on  the  phugoid  motions.  While  each  artificial  derivative  reduces  the  period 


to  62.8  sec. 


c 


M, 


decreases 


somewhat,  and  £ 


causes  a 


Pon 

slight  increase  over  the  normal  damping.  The  pitch  angle  response  per  unit 
is  raised  to  1.81  by  and  lowered  to  0.90  by 


'A Dll 


The  phugoid  damping  can  be  increased  to  J?6.1#  critical  by  addition  of 
Cw...  > or  ■ as  seen  Figure  (11),  All  three  artificial 


C 

4 * '—'‘rr q 

derivatives  maintain  the  period  close  to  the  normal  88,6  sec.  value,  A significant 
increase  in  angle  of  attack  response  is  noted  with  the  use  of  or  Cj# 

The  pitch  angle  lag  is  reduced  to  5h°  with  C» t and  £2°  with 

increases  the  lag  to  128°, 


Du 

while 


e 


The  usefulness  of  these  phase  diagrams  is  apparent  from  the  ease  with 
which  the  amplitude  and  phasings  of  the  motions  can  be  obtained.  It  should  be 
noted,  that  these  diagrams  can  be  obtained  quite  readily  by  graphical  analysis. 

Besides  enabling  the  theoretical  analyst  to  gain  a physical  insight  into  the 
motions,  these  diagrams  can  be  used  to  predict  the  transient  response  of  other 
variables,  once  the  response  of  a single  variable  has  been  determined.  Thus, 
assume  that  the  D © response  to  a particular  input  has  been  determined  as: 

D Q = c, 

part. 

where  D0  part,®  particular  solution  for  'D©  • The  response  to  any  other 
variable  X can  be  determined  as: 

X «=  Axj  c,  e ^ a,os(pi  tr-f  <Jy  -f 

-j~  X part. 
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where 


f\x  1 ■ 

fixi  “ 


amplitude  ratio  of  X to  QQ  in  short  period 
( c short  period  frequency) . 

amplitude  ratio  of  ^ to  p 0 in  phugoid 
b phugoid  frequency) 

phase  lead  of  X from  Qg  in  short  period 
phase  lead  of  X from  q Q in  phugoid 
particular  solution  for  X, 


The  new  particular  response  to  the  input  disturbance  must  bo  calculated 
separately.  In  the  case  of  step  disturbances  or  initial  conditions  this 
particular  response  can  be  found  for  combined  phugoid  and  short  period  analysis 
from  the  formulae: 


Step  Elevator  Particular  Solutions 


^ part. 

^part. 

“"’'part. 


See  formulae 
state  values 

^ss 

in  Appendix 

for  steady 
on  Page  43 


~Jr  part 


O 


Particular  Solutions  When  Input  is  Initial  Condition  on  u or  c*r 


part. 


6. 


part. 


°Vrt.C  * 


£ P^t, 


0 


With  the  particular  solutions  known,  the  response  in  the  desired 
motion,  X j can  be  readily  determined  with  the  amplitude  and  phase  ratios 
(Ax  and  x ) determined  at  both  short  period  and  phugoid  frequencies. 
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RESULTS  AND  CONCLUSIONS 

jr 

If  On  the  basis  of  this  theoretical  investigation  to  provide  a wide 

range  of  artificial  stability  and  control  to  the  F-9U  airplane,  the  following 
. it  results  and  conclusions  can  be  noted: 

I ) o 

1,  Use  of  elevator  signals  proportional  to  o?  and  will 
provide  a wide  range. of  short  period  frequency  and  damping 
control. 

2«  Usp  of  auxiliary  surface  position  signals  proportional  to 
i ; /I  q/q  and  u/q  will  provide  a wide  range  of  phugold  period 

( ( and  damping  control. 


II 


I 


i 


3*  Use  of  elevator  position  signals  varying  with  stick  position 
will  change  the  fixed  stick  stability. 

ij.  Use  of  stick  force  signals  varying  with  stick  position  and 
angle  of  attack  will  change  the  free  stick  stability  of  the 
airplane . 

9.  Seven  artificial  stability  derivatives  will  be  provided  on 
the  F-9U  with  a hydraulic  servo  system  on  the  elevator 
and  stick,  and  with  an  electric  servo  motor  on  the  auxiliary 
surface, 

6,  An  auxiliary  surface  of  0.U3  sq.  ft.  will  be  constructed  and 

installed  in  the  nose  section  for  phugoid  control.  The 
surface  will  be  geared  down  to  2 deg/ sec.  to  avoid  normal 
acceleration  peaks  due  to  forward  velocity  gusts  with  w 
present.  ^ 

7,  The  phasings  of  the  phugold  and  short  period  motions  can 
be  investigated  readily  by  vector  phase  diagrams. 

8,  The  phasings  of  the  oscillatory  modes  may  be  an  important 
parameter  in  the  pilot's  evaluation  of  a flight  configuration. 
Further  insight  into  the  effect  on  the  human  pilot  of  the 
phasings  in  the  short  period  and  phugoid  response  should  be 
gained  by  analysis  of  flight  traces  as  well  as  by  theoretical 
studies. 


i 
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TABLE  I 


F-9UA  NORMAL  STABILITY  DERIVATIVES 


Constants: 


e * 6.72  ft. 

.525 

JH_ 

- - 0.90 

K 

S - 237.6  ft. 2 

d c* 

1+7.83  ft. 2 

q 

« .720 

ft. 

St  - 

ce 

**  1.03  rad./ft. 

*e  * 

8.70  ft. 2 

Variables: 

Indicated  Airspeed,  (mph) 

135 

219 

292 

• 290 

365 

Mach  Number 

.18 

.1+3 

.57' 

.70 

.71 

True  Airspeed,  V (mph) 

135 

300 

l+oo 

1+71+ 

500 

Altitude  (ft.) 

S.L< 

, 20.000 

20,000 

30,000 

20s000 

C.Q.  (gMAC) 

28 

29 

29 

29 

29 

-.01+9 

-.055 

-.055 

-.055 

-.055 

G.W.  (lbs.) 

12,359 

13,611+ 

13,611+ 

11,365 

13.611+ 

I (slug  ft.2) 

26,51+5 

26,51+3 

26,51+3 

22,158 

26,51+3 

iB 

6.12 

5.51+ 

5.51+ 

5.51+ 

5.51+ 

2 ijc 

5.23 

5.21 

5.21 

5.21 

5.21 

S f (deg.) 

1+5 

0 

0 

0 

0 

CL 

1.12 

.1+66 

.261+ 

.222 

.168 

(1/rad) 

5.12 

5.51 

5.90 

6.1+8 

6.hl 

CD 

.135 

.0307 

.0235 

.021+2 

.0237 

CDp,  (1/rad) 

.653 

.261+ 

.11+7 

.11+0 

.110 

Gm^  (1/rad) 

-.251 

-o30U 

-.321+ 

-.357 

-.353 

Cm  n (1/ deg) 

a* 

-.0155 

-.0165 

-.0165  - 

-.0165 

-.0165 

(1/rad) 

-8.96 

-9.28 

-9.28 

-9.28 

-9.28 

%«  (1/rad) 

— 1+  .10 

-1+.21+ 

-1+.21+ 

-1+.21+ 

-1+.21+ 

'Z'  (sec.) 

3.1+5 

3.20 

2.1+0 

2.1+1 

1.92 

203 

1+19 

1+19 

1+97 

kl9 

ch  s ft/rad) 

-.1+30 

-.1+78 

-.1+82 

-.1+60 

-.1+55 

0h  (l/rad) 

-.15 

-.10 

-.09 

-.09 

-.06 
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TABLE  V 

CONTROL  STICK  PARAMETERS  OF  THE  NORMAL  F-9 


(deg/g) 


Push  Downs  Turns 


(lbs/ g) 


Push  Downs  Turns 


Vi  = V* 


Vi  -=  V*. 


(deg/mph  ind)  (lbs/mph  ind) 


-.01*9 

135 

lu36 

U.36  ? 
+.822/n 

U.35 

U.36 

+1.39/n 

.0521* 

.01*38 

-.055 

219 

1.879 

1.879 
+ . 3Il*/n2 

8.06 

8.06  „ 
+1.66/n2 

.011*1* 

.0585 

-.01 

292 

.338 

.338 

+„l?7/n2 

.977 

0.977  0 
+1.6 9/n2 

.00109 

-.001*88 

-.095 

292 

1.059 

1.059 

4.177/n2 

8.60 

8.60 

+l,69/n2 

.00601 

.01*71* 

-.10 

292 

1.782 

1.782 

+.177/n2 

16.21* 

16.21* 

+1.69/n2 

.0109 

.0996 

-.055 

290 

.923 

.923,, 
+ ,179/n2 

7.21 

7 = 21  o 
+l»6l/n2 

.00512 

o0386 

-.055 

365 

.677 

0677 

4 .113/n2 

9.01 

9.01 

+1.65/n2 

.00309 

.01*03 
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TRANSIENT  RESPONSE  TO  STEP  ELEVATOR 
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TRANSIENT  RESPONSE  TO  STEP  ELEVATOR 
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FIGURE  3 

TRANSIENT  RESPONSE  TO  STEP  AUXILIARY  SURB&CE  DEFLECTION 

PHUGOID  EFFECTS  ,Cmj  -02 9 ,CmDl/ia  .171 


n*vg*s 
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NORMAL  ACCELERATION 


aV  ~ mph 


FORWARD  VELOCITY 
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TRANSIENT  RESPONSE  TO  FORWARD  VELOCITY  GUST 

NORMAL  AIRPLANE 


NORMAL  ACCELERATION 


n ivj  9 


DEGREES 


ANCLE  OF  ATTACK 


DEG. / SEC. 


PITCH  RATE 


& 

DEGREES 


PITCH  ANGLE 


A V-v 

mph  Ind. 


FORWARD  VELOCITY 


-£  -v.  SECONDS 


CONDITIONS 
yL  = 292  mph 

H - to, ooo  ft. 

/d  CL  - “*05i 

O,  W.  - 13, 614  Ibi 


WADS  52-P)t8 


F-94  FIGURE  5 

TRANSIENT  RESPONSE  TO  FORWARD  VELOCITY  OUST 

Cm." .029 
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SHORT  PERIOD  PHASING 
FREQUENCY  VARIED 


FIGURE  8 
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SHORT  PERIOD  PHASING 
DAMPING  VARIED 

90° 


normal  f*m  concitionb 


APPENDS 


VECTOR  ANALYSIS  OF  MOTION 

A force  vector  diagram  can  be  drawn  for  any  equation  when  the  direction 
and  amplitude  of  the  variables  is  known  at  any  natural  frequency  of  the  motion. 

In  tha  aircraft  longitudinal  motion,  the  phugoid  and  short  period  modes  are  of 
greatest  interest.  The  amplitude  and  phase  relationships  of  the  variables  will 
be  determined  at  these  frequencies. 


Short  Period 


The  short  period  mod®  can  be  represented  by  the  equations} 

(d  -f  • — OB  — 0 


*5f 


l 


* Hr 


j>  - 


0 = 0 


The  natural  undamped  frequency  and  damping  of  the  short  period  can  be  obtained 
from  Equation  (1)  and  froa  the  relation: 

t*^>L  — 


the  relation  between  Os*  and  can  be  expressed  as: 

Do*  ■ /so0  * 


where  ^ 

€ 0 sm 

This  expression  can  be  found  in  Reference  (3). 

The  lift  equation  can  now  be  used  to  obtain  D © in  terms  of  cK  ‘ 

°<r 


A graphical  solution  will  yield  sufficiently  accurate  results.  The  presence  of 
Oh  in  the  short  period  can  be  predicted  accurately  from  the  drag  equation: 

( D 4-  C0  ) « 4-  ( ^r*~ 

Since  ^-0  is  very  small  compared  to 
this  equation  becomes: 

a<~  - (c-^~  ~c.L)^ 

which  can  be  solved  graphically. 

The  quantities  D © and  0 U have 
relation  to  unity  at  zero  phase.  These 
below  for  the  normal  airplane. 


— CL)<*  -f- 

at  the  short  period  frequency. 


been  determined  as  vectors  in 
quantities  are  shown  graphically 


PQ 


id* n / 9Q°  * £d  + 


Cl 


S7  _ 


SHORT  PERIOD  PHASING 
FREQUENCY  VARIED 


30* 


phugoid 

The  phugoid  mode  can  be  represented  by  the  three  equations  of  motions 
{ 0 T'vpj  to  - Cl~)  * +■  0=0 

uu  +-  ( £)  -t-  ) < - o & *=  0 

-+  [d2- 6 = o 
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If  0^6  Is  neglected,  the  latter  two  equations  can  be  combined  as: 


-4-Cv  -S&dpt*  *+  —JLo -Sia^^2s.o< 


The  vector  Qc^  . k/w  ay  J 90°  +£i 


(ip'*)* +[:£<. 


‘f'"^ 


311(1  - lJh  ^ /a DS-igp 

The  natural  undamped  phugoid  frequency  and  damping  can  be  obtained  from  Equation  (U) 
and  from* 


''  PtrJPpi 

e o ■ snrl 


a*i  can  then  be  determined  in  terms  of  U,  . The  drag  equation  can  then  be  used 
to  determine  graphically  the  vector  properties  of  © in  terms  of  unit  at 
aero  phase.  The  normal  F-9h  phugoid  motions  are  shown  vectcrially  below: 


PHASING  OF  THE  FHUGOID 
PERIOD  VARIED 


O 
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ANALOG  COMPUTER  SETUP 


The  longitudinal  equations  of  motion  for  the  F~94  airplane  at  I4OO  mph, 
20,000  ft.  altitude  for  $»$%  static  margin  ar6  listed  below 


du 

dV«i 


,049  u +.122**  -,,276  0 


- * -.56  u -6»14®f  + 


4.M,  ^,2-  -107“>  -1-60  M -3-6^  'su  Se- 

where  fcj  - 5.0  sec. 

The  elevator  was  positioned  by  the  signals  1 
S.  — S * .0219  -.0307  u -0876^, 


, = * -0S19  % -°307  U -0876% 
applying  .AAh  - -.389,  c -.01*3,  Cy*  * .029, 

Oar  ” i 


Inputs 

A forward  velocity  gust  is  applied  to  the  airplane  on  the  analog  by 
entering  an  initial  condition  on  CL  . This  is  equivalent  to  a forward  velocity 
increment  of  the  air  relative  to  the  ground.  Since  the  plane's  velocity  with 
respect  to  the  ground  will  not  change,  this  gust  will  give  a sudden  increase  in 
the  velocity  of  the  plane  with  respect  to  the  air  as  seen  mathematically  below: 


v - L{ 


pg  + uag 


This  aag  provided  the  disturbance  to  excite  the  phugoid  oscillation. 

It  should  be  noted  that  a sudden  increase  in  *^pa  will  cause  an 
impulse  in  the  rate  of  change  of  airspeed  as  measured  by  a pitot  tube.  Thus, 
as  is  provided  by  sensing  the  rate  of  change  of  ^ ' pa,  an  impulse 

in  pitching  moment  would  be  applied  to  the  airplane,,  The  auxiliary  surface  is 
limited  to  2 deg/sec.,  however,  so  that  the  rate  of  change  of  pitching  moment 
is  drastically  limited. 

In  order  to  examine  the  response  of  the  airplane  without  this  velocity 
limiter,  initial  conditions  on  U and  0 were  entered  on  the  analog.  The  condition 
on  0 arises  from  the  previously  noted  effect  o£Cyyi(?u  . Thus,  the  integral  of 
the  moment  equation  must  be  in  balance  before  the  gust  is  applied  at  time  -0  and 
just  after  application  at  time  +0,  This  yields, 
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as  40  ca 


This  relation  between  and  U,  is  necessary  to  represent  u gust 
properly  when  Oy*  of  #1TL  is  present,  <{■ 


An  angle  of  attack  gust  or  a vertical  velocity  gust  will  be 
represented  similarly  by  initial  conditions  on  g*  and  0 , Thus,  for  the 
normal  airplanes 


£,  G 


-1*60 


When  AC-w,  . , 

101  Ocx  a 


is  present,  an  of,  gust  creates  a larger  as: 


if  AC>np  ■ =-0U3j  ( &/<ktMo  = -8. bo  o<Q 


0 « 

Besides  the  inputs  of  initial  conditions  on  U , ,0  and  e>< , 0 ; 

the  airplane  response  to  a step  elevator  and  to  a ramp  elevator  was  examined. 
The  response  to  a ramp  represented  the  most  severe  response  possible  with  the 
velocity  limiter  used  on  the  auxiliary  surface.  It  should  be  noted  that 
although  the  equations  are  s^t  up  for  elevator  inputs,  auxiliary  surface  inputs 
ean  bo  obtained  by  ratioing  the  response  in  terms  of  C Irntr  /Cw- 

Wiring  Diagram 

The  wiring  diagram  is  included  next  to  indicate  the  circuitry 
necessary  to  obtain  the  responses* 
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STEADY  STATE  VALUES 


The  following  formulae  are  listed  to  aid  in  determining  the  steady- 
state  values  of  the  aircraft  motions  following  a step  elevator  deflection: 


CK  - — 

CUef 


®r  - w 


L.  ' 


Ltcd&±,k lL 

l fi6  £•  j.  / (assuming  <~s_s  * 0) 

The  steady  state  values  in  the  short  period  mode  of  motion  are: 


e 


C, 


* - c. 


1 


i 2/^/ 


o'c1 


ft. 


2*tf 


/ 


The  maximum  value  of  «=><  in  the  short  period,  and  of  ^ ^ and  G in  the 
phugoid  can  he  computed  from: 

max. » ^ steady  state  (!•».  % overshoot  ) 
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where  either  oj  N^u'j  or  © can  be  represented  by 
The  percentage  oversnoot  is  a function  of  ^ along  and  is  sh 


shown  below: 
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